Current Fluctuations and Electron-Electron Interactions in Coherent Conductors 
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We analyze current fluctuations in mesoscopic coherent conductors in the presence of electron- 
electron interactions. In a wide range of parameters we obtain explicit universal dependencies of 
the current noise on temperature, voltage and frequency. We demonstrate that Coulomb interaction 
decreases the Nyquist noise. In this case the interaction correction to the noise spectrum is governed 
by the combination T n (T n — 1), where T n is the transmission of the n-th conducting mode. The 
effect of electron-electron interactions on the shot noise is more complicated. At sufficiently large 
voltages we recover two different interaction corrections entering with opposite signs. The net result 
is proportional to ~}2 n T n (T n — 1)(1 — 2T„), i.e. Coulomb interaction decreases the shot noise at low 
transmissions and increases it at high transmissions. 
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I. INTRODUCTION 

Recent advances in nanotechnology enable detailed in- 
vestigations of a variety of quantum effects in mesoscopic 
conductors. These investigations are of primary inter- 
est because of fundamental importance of such effects 
as well due to rapidly growing number of their poten- 
tial applications. A great deal of information is usu- 
ally obtained from studying electron transport. Addi- 
tional/complementary information can be extracted from 
investigations of fluctuation effects, rfijwh as shot noiseE 
For instance, it was demonstratedau'cl that the power 
spectrum of the shot noise in coherent mesoscopic con- 
ductors is expressed in terms of the parameter 



(1) 



Here and below T n stands for the transmission of the n- 
th conducting channel of a coherent conductor. Thus, 
since transport measurements only allow to determine 
the combination 



1 
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(2) 



studies of the shot noise provide additional valuable infor- 
mation about the transmission distribution of conducting 
modes. 

The above results apply to the situations when inter- 
action between electrons can be neglected. In the pres- 
ence of electron-electron interactions the Landauer con- 
ductance (|J) and the I — V curve are modified in a non- 
trivial way. Recently it was shownQ that the I —V curve 
of a (comparatively short) coherent conductor with ar- 
bitrary transmission distribution T n in the presence of 
interactions can be expressed in the form 



R^L = l-f3f(V,T), 



(3) 



where f(V,T) is a universal function to be defined be- 
low. This result holds in the limit of large conductances 



R <C R q — h/e 2 or, otherwise, at sufficiently high tem- 
peratures/voltages. It demonstrates that the magnitude 
of the interaction correction is controlled by the same 
parameter (3 (|l|) which is already well known in the the- 
ory of shot noise. Physically this result can easily be 
understood since both phenomena are related to discrete 
nature of the electron charge. Hence, there exists a direct 
link between shot noise and interaction effects in meso- 
scopic conductors. 

It is obvious that not only the I — V curve (Q) but 
also shot noise as well as higher moments of the current 
operator should be affected by electron-electron interac- 
tions. This paper is devoted to a detailed investigation 
of current fluctuations in mesoscopic coherent conductors 
in the presence of electron-electron interactions. Previ- 
ously various aspects of this problem have been studied 
for a particular case of tunnel junctions in the Coulomb 
blockade regime, see, e.g. Refs. ||7|J||| The effect of 
interactions on shot noise in 2d diffusive conductors at 
sufficiently high temperatures was recently addressed in 
Ref. 

S2Efi_we will employ a model of a coherent conduc- 
Ej£3. Within this model we will demonstrate that 
interactions lead to two different corrections to the shot 
noise spectrum. One of these corrections scales with the 
parameter (3 (|l]) . This correction is always negative, sim- 
ilarly to that found in Eq. (||) for the I — V curve. It 
describes (partial) suppression of the current noise due to 
Coulomb blockade. In addition to this correction we shall 
find another one, which is proportional to the parameter 



7 : 



12n Tn 



(4) 



This second correction is positive, i.e. it leads to rela- 
tive enhancement of the shot noise. The latter correction 
turns out to be important only at voltages exceeding both 
frequency and temperature and is negligible otherwise. 
Thus, at sufficiently high voltages two interaction correc- 
tions - negative and positive - compete, for 7 > (3/2 the 
second one wins and, hence, in this case an overall en- 
hancement of the shot noise by interactions is predicted. 
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FIG. 1: The circuit under consideration. The scatterer (de- 
noted by a cross) has a capacitance C and is connected to the 
voltage source V x via an impedance Zs- 



Our paper is organized as follows. In section 2 we 
will describe the model and highlight our key results. A 
detailed derivation of these results is then performed in 
section 3. Our main conclusions are briefly summarized 
in section 4. Most of the technical details, such as the 
derivation of our effective action as well as a few other 
issues are presented in Appendices A,B and C. 



II. THE MODEL AND KEY RESULTS 

Similarly to Rcfs. ||,[ll],[l2] we will consider a coherent 
scatterer between two big reservoirs. The scatterer is 
described by an arbitrary distribution of transmissions T n 
of its conducting modes, and the corresponding Landauer 
conductance 1/R is defined in Eq. (||) . The transmissions 
T n are assumed to be energy independent. Phase and 
energy relaxation may take place in these reservoirs but 
not inside the scatterer, i.e. the scatterer is assumed 
to be shorter than both dephasing and inelastic lengths. 
The scatterer region has an effective capacitance C. For 
simplicity we will also assume that the charging energy 
Eq = e 2 /2C does not exceed the typical inverse traversal 
time. The scatterer is connected to the voltage source V x 
via linear external impedance Z$ (see Fig. 1). Here 
we restrict ourselves to a simple case Zs{u>) = Rs- If 
necessary, generalization of our calculation to arbitrary 
Zs(uj) can be performed in a straightforward manner. 

In what follows we will investigate the current noise 
and evaluate the correlation function 

s(t,t') = \{i{t)i{t') + i{t')i{t)) - (i) 2 , ( 5 ) 

where / is the current operator in the circuit of Fig. 1. 
This correlator can be expressed in the form 



S(t,t') =S ni (t,t') +SS(t,t'), 



(6) 



where S 
rent nois< 



he noninteracting contribution to the cur- 
and 5S is the correction due to electron- 



electron interactions inside the scatterer. This correction 
will be evaluated in the most interesting "metallic" limit 



.9o = 9 + 9s > 1. 



(7) 



Here we introduced dimensionless conductances of the 
scatterer and the shunt, respectively g — R q /R and gs = 
Rq/Rs- Eq. (Q) implies that at least one of these two 
dimensionless conductances is required to be much larger 
than unity. 

Quite obviously, the latter correlator (||) should depend 
on both R and Rs- We also introduce another correla- 
tor S(t, t') defined by the same Eq. (||) in which one 
should substitute the current operator across the scat- 
terer / — > I sc . The two correlators S(t, t') and S(t, t') are 
not independent. With the aid of the current conserva- 
tion condition and performing the Fourier transformation 
with respect to t — t', one easily finds the relation 



|f(l + ^i? 2 C 2 )5, 

R S , ,2 D 2^2\ • - 



R- 



(l+uo z R z C z )ujcoth. 



2T' 



(8) 



where i?o = RRs/(R + Rs)- The second term is due the 
noise produced by the external resistor Rs which has to 
be subtracted in order to arrive at S u . 

In general also the correlator S u depends on both R 
and Rs- However, in the limit Rs ^> R the dependence 
on the shunt resistance is weak and can be neglected. In 
this case the interaction correction to the current noise 
spectrum SS^ depends only on the properties of the scat- 
terer. Below we will present our key results for SS^ only 
in this limit. More general expressions can be found in 
Sec. 3. 

Let us define the average voltage across the scatterer 
V = V X R/(R + Rs) and consider first the limit of rela- 
tively small voltages. At sufficiently large temperatures 
and/or frequencies we find 

*& = -^£, ifT>^ c ,|eV|,M, (9) 



3R 

I3E C 
R 



, j£\u\x>T,gE c ,\eV\. (10) 



At lower temperatures and frequencies we obtain 
<J& = -^ln^, ifM,|ey|«T« 5 £ c , (11) 



Ra 



T 



5S u = -^M\n^, ifT,| e y|«M« 3 £ c . (12) 

Rg M 

These results apply as long as either temperature 
or frequency exceeds an exponentially small parame- 
ter gEc exp(j=jg/2). For even smaller frequencies and 
temperaturestil we get 



5S^ 



0u> u> 
li COth 2T- 



(13) 
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Note that the above expressions could also be antici- 
pated from the fluctuation-dissipation theorem (FDT) 
combined with the resultsEl. Indeed, in the limit of low 
voltages the current noise is described by the standard 
Nyquist formula. Hence, in order to satisfy FDT one 
should simply substitute the effective conductance (||) 
into this formula. In this way one gets the interac- 
tion correction SS U proportional tp f3f. For instance, in 
the low frequency limit one findsB f(0,T) ~ Ec/ST for 
T » gE c , f(0,T) ~ (2/g)hx(gE c /T) for exp(- 5 /2) « 
T/gEc < 1 and f(0,T) ~ 1 for T < gE G exp(-g/2). 
Combining these expressions with FDT one immediately 
reproduces Eqs. (||), ( fiT| ) and (13). 

It is worth stressing that here we evaluate the current- 
current, correlation functions directly and do not use the 
resultsa together with FDT. However, it is satisfactory 
to observe that FDT is explicitly maintained in our cal- 
culation and the results derived here are fully consistent 
with those of Ref. || 

Now let us turn to the case of relatively large voltages 
V where the shot noise becomes important. As it was 
already announced, in this case the correction to the noise 
power spectrum is proportional to the parameter 



(3-2 1 



£ n T„(l-T„)(l-2T n ) 



In particular we obtain 



6S,„ = - 



((3 - 2j)E c 
R 



R q \eV\ ' 

ifT.M « \eV\ <ZgE c , 

ti\eV\Z>T,gE c ,\w\' 



(14) 

(15) 
(16) 



We note that this correction can be either negative or 
positive depending on the relation between the parame- 
ters P and 7. Thus, in contrast to the limit of low voltages 
(Nyquist noise), one cannot conclude that shot noise is al- 
ways reduced by interactions. This reduction occurs only 
for conductors with relatively low transmissions (3 > 27, 
while for systems with higher transmissions the net ef- 
fect of the electron-electron interaction enhances the shot 
noise. In the important case of diffusive conductors one 
has (3 = 1/3, 7 = 2/15 and, hence, 



1 
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In this case the shot noise is reduced by interactions. 

The above results have a transparent physical inter- 
pretation. At low voltages the power spectrum of the 
Nyquist noise is proportional to the system conductance 
tx J2 n Tn- Since in the presence of interactions the con- 
ductance acquires a correction proportional to j3, the 
interaction correction to the Nyquist noise should scale 
with the same parameter (Q). On the other hand, shot 
noise is determined by the combination $^ n (T n — T%). 
Accordingly, the interaction correction to the shot noise 



power should consist of two contributions. One of them 
comes from E n T„ and is again proportional to (3. An- 
other contribution originates from the interaction correc- 
tion to E ra which turns out to scale as 27. Since these 
two corrections enter with the op pos ite signs we imme- 
diately arrive at the combination (J14|). 

We also point out that the third cumulant of the-cur- 
rent operator for noninteracting electrons is knowr£3 to 
be proportional to the parameters ([!]) and ( |l4| ) respec- 
tively at low and high voltages. Following the same ar- 
guments as above we can conjecture that the interaction 
correction to the third cumulant should scale as (3 — 27 
at low voltages, while in the limit of large voltages one 
can expect that this correction is governed by the com- 
bination [3 — 67 + 65, where 



En^n(l-r n ) 
En Tn 



(17) 



This conjecture can also be generalized to higher cumu- 
lants of the current operator. 

We would like to emphasize that - although the above 
conjecture seems intuitively appealing - it should still be 
verified by means of a rigorous calculation which is be- 
yond the scope of the present paper. In the next section 
we will concentrate on the current noise and will provide 
a detailed derivation of the results presented above. 



III. EFFECTIVE ACTION AND CURRENT 
NOISE 

Similarly to Ref. [| we will use the effective action tech- 
nique in order to evaluate the current-current correlator 
for the system depicted in Fig. 1. It is convenient to 
introduce the quantum phase variable ip which is pro- 
portional to the integral of the fluctuating voltage (see 
Appendix A). We will proceed within the Keldysh for- 
malism and introduce two phase variables tpi^ related 
to the two branches of the Keldysh contour. Defining 
ip + — (ipi + v?2)/2 and tp~ — fi — if2 one can denote the 
overall phase jumps across the scatterer as tp + + eVt and 
ip~ . Correspondingly, the phase jumps across the Ohmic 
shunt are (eV x — eV)t — ip + and — tp~ . The symmetric 
current-current correlation function (||) can be expressed 
as follows 



\{i(t)i{t') + i(t')i{t)) 



+ 



6<pg(t)6<ps(t') ' 4S<p+(t)6<p+(t') 



Vi P ± 
Stod^ ( 18 ) 



see Appendix A for further discussion. By we denote 
the phase jumps over the ohmic shunt. The variational 
derivatives in Eq.(|l^) act on the shunt part of action. 
Here 5 to t[ ( / : ' ± ] is the total action of our system 



Su A y k \ = s[<p k ] + s 8 y k ], 



(19) 
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where the term 
iS s [f ±] - 



e 2 Rs Jo 



dttp-(t) 



1 



, , dh / dt 2 a(h - t 2 )ip (h)p (t 2 ) 
2e Rs Jo Jo 

comes from the shunt, S^^] is the scatterer action and 

2 



a(t - 1') 



1 



ttT 



tt Vsinh[7rT(t - t')] 



(21) 



A detailed derivation of the action S^y^] is carried out 
in Appendices A, B and C. The main idea of this deriva- 
tion is to expand S[<p ] in powers of ip~ keeping the full 



nonlinearity of the corresponding terms in (p + . This pro- 
cedure is just the quasiclassical approximation for the 
phase variable. It is parametrically justified under the 
condition ((?]). In Ref. || the action /Sfy^] was evaluated 
up to the second order in (p~. This is sufficient to de- 
rive the current-voltage characteristics of the scatterer. 
However, in order to describe the current noise it is nec- 
essary to expand the action <S , [(/? ± ] further and to retain 
all terms up to the third order in ip~ 



(22) 



This expansion is analyzed in Appendix C. We will now 
use these results and explicitly evaluate the current- 
current correlator (ph. 



A. Contribution of first and second order terms 



Let us first restrict our attention to the contribution of the first and second order terms in (E2). They read 



e Jo 



oo poo 



dh / dt 2 a{h - t 2 )<p- (h)cp-(t 2 ) {1- (3 + cos [eVih-^ + tp+ih) -<p + (t 2 )]} 
le Jo Jo 



Employing Eq. (|18| ) we obtain 



\{I{t)i{t') + I(t')i{t)) = + e 2 (K(t)K(t<)), 



where we defined 



K(t) 



e 2 Rs 



eV x - eV - ip + + i / dta(t - t)cp" (t) 
o 



(23) 



(24) 



Note that here we neglected terms which originate from the second variational derivative with respect to ip + in ( |l8| ) 
since these terms are smaller in the parameter than those kept. Angular brackets in Eq. ( p3| ) imply averaging 
with the path integral 



J Vi P ± {...) exp (iS<U [^] + iS^l^} + iSsl^} 



Rewriting the correlator (KK) as 



W (0) = -Jn/V^ i)j(t 



-expjz ^ S^l^+iSslv^+iJ^ 



dir)(t)K(i) 



(25) 



and performing a shift of (p u — ► (p^ + r) u Ro/ R S {1 — iloRoC), we obtain the expression for the Fourier transformed 
noise spectrum (|^) 



w I 1 

Su, = uj coth — < Re ■ 



(3 \ P[a(t-t')(cos(eV(t-t') + ^+(t)-^+(t')))} 



2T Z(lj) Rntt 



Rntt 



(26) 



-^[ m r-(H(t)H(t'))] i 



where 



Hit) = I dta(t - i)cp-(i) [l - cos (eV(t - i) + cp + (t) - ip + {t))] . 
'o 



(27) 
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In Eq. (t26J) we also introduced the following notations 



Z(uj) 



Rs 



1 



R- 1 - ituC 



n 



(R + Rs) 2 
R 2 



(28) 



Making use of the relation (^) , we arrive at the correlator S u , 

S u = i ((1 - /3) W coth ^ + (3[a(t - t')(cos(eV(t - t') + <p+{t) - v + {t')))]S) + ■ ■ 



(29) 



where . . . stands for the terms containing [(7J(i)) 2 — (H(t)H(t'))] in Eq. (26). 

Within our approach interaction effects are described by the terms containing the fluctuating vari 
formally sets this variable equal to zero, from (Eq) one immediately recovers the noninteracting resul 



If one 



3? = (!-/*)£ -th £ + ± eV) coth 



(30) 



Taking the phase fluctuations into account we arrive at the expression for the interaction correction to (|30j). However, 
the corresponding expression turns out to be incomplete in two respects. First, one of the terms does not satisfy 
FDT. Second, the correction to (^) obtained in this way scales with the parameter (5 in both limits of small and large 
voltages. While in the former limit (Nyquist noise) this result is understandable and consistent with Ref. [|, a t large 
voltages (shot noise) one also expects an extra contribution. Its existence can be anticipated because the shot noise 
is governed by the combination ^ n T„(l — T n ) and not simply by ^ n T n as the Nyquist noise, see also our discussion 
in Sec. 2. 

Both these problems are remedied by taking into account the third order in tp~ contribution to the effective action. 
This will be demonstrated in the next subsection. 



B. Corrections due to third order terms 



Following the analysis in Appendix C we identify two different contributions to the third order term 



The first contribution has the form 



tb P [V J _ 6e*R 



dr for(r))V(r) 



(31) 



Taking this term into account and repeating the above analysis we arrive at an extra contribution to the current 
noise in the form —e 2 (35(t — t')/2RC, see also Appendix C. Adding this contribution to Eq. (26) and subtracting the 
noninteracting result (Bfj) we arrive at the interaction correction 



6SW = 



RkQ, 



t(t - t') cos (eV(t - t')) [e- F{t - t ' ) - l) - e 2 S(t - t')/2C 
The function F{t) results from averaging over the phase fluctuations 

{cos[eV(t - t') + tp + (t) - tp + (t')}) = cos[eV(t - t')}er F ^ t - t '\ 

This function has the form 



F(t) = e 2 R 2 



dcu 1 — cos tot 
2^ ljHI 



— ) u> coth h ~S^(uj ± eV) coth 







Ro R 



uj (3 



2T 2R 



u±eV 
2T 



(32) 



(33) 



(34) 



We also note that in Eq. (|32j) we omitted the last term of Eq. (|26j) which contains averages of the function H (|2 
Our analysis demonstrates that these terms are small in all the regimes considered below. 
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What remains is to evaluate the correction to the shot noise from the second contribution to S^ 3 \ The derivation 
of this contribution is presented in Appendix C. Here we only quote the result: 



ni'yT 
6e 2 R 



dyi 



dy3 



v> (yi)f (2/2)95 (2/3) 



sinh[7rTj/2i] sinn^T^] sinh[7rT?/i3] 



{sin (eVy 2 i + </? + (j/2) - <P + (yi)) + sin (eVy 32 + V + (y3) - t f + {y2)) + sin (eVyi 3 + ip + (yi) - <p + (y3))} 



(35) 



where yij = yi — yj and the parameter 7 is defined in Eq. (Q). 

At the first glance this contribution to the effective action could be considered unimportant. This is indeed the case 
in several limits. For instance, at sufficiently small transmissions j3 S> 7 the term (pq) can obviously be neglected. 



In the limit of low voltages one can, making use of the condition (]?]), expand Sj in small phase fluctuations tp^. 
Then one gets Sj proportional to the combination ((,2 + ) 3 (<y9~) 3 which can be dropped as compared to other terms 
provided go 3> 1. However, in the limit of large voltages the term (|3^) gains importance and - as we shall see - 
provides significant contribution to <56> w . 

Proceeding along the lines with the above analysis we find 



2tt 2 7T 3 
dt 



tlx 



(f - e -*ARoC) ( C oa[eVt] - cos[eVx}) cos ut 



sinh[7rTx] sinh[7r7Y] 
This expression will be analyzed below in Sec. 3F. 



f 



f 



sinh[TrT(a; - t)] sinh[7rT(x + t)} 



r 



(36) 



C. Relation to FDT 

Before we proceed with the analysis of the above re- 
sults let us establish some useful general expressions and 
illustrate the relation between our approach and FDT. 
Throughout this subsection we will only consider the 
limit of small voltages eV <C 1/ RqC and neglect the de- 
pendence of the function F(t) on V. In the spirit of the 
P(£ , )-theoryt3 let us define the function 



P(E) 



dte iEt e -*{t) } (37) 



#(t) = F(t)\ v=0 + 



ie 2 R Q 



sign[t] ( 1 



-\t\/R C 



This function obeys the "detailed balance" symmetry 
P(—E) = e~ E / T P{E) which follows from the property 
<f>(i— (i/T)) = $(— t). Let us also introduce the function 



A/L, 



1 



AttRq 



E 



± 



dE 



E(l 



-(uj±e 



1 - e- E / T 
and rewrite it in the form 



P(uj±eV-E), (38) 



1 

Ro J-c 



a(t)e~ F W cos(eI//)cos 



dte 1 



-0 

9o v 



2C 

-t/RoC 



S(t) 



(39) 



We observe that, since in the interesting for us limit g 
1 the argument of cos 



the accuracy ~ 1/So one can use the function ( |39| ) in 
order to analyze the result (|32"|). 

Proceeding further let us rewrite Eq. (pq) as 



AAH = i^coth^^X( W ±ey) 



where 



1 



-u/T 



dEE 



(40) 



tP{u-E). (41) 



After a simple algebra from Eq. ([H|) we obtain I(i 
jfc + 81{uj) and 

2 

Ro 

-t/RoC 



dtsm(tot)e F ^a(t)sm 



90 



1 



(42) 



Comparing the above expressions with Eq. (|32| ) we arrive 
at the following correction to the current noise 



0*° Y J C0th^^8I{u±eV). (43) 



2RnVl 



2T 



In order to illustrate the relation between our results and 
FDT we notice that in the relevant limit go 3> 1 the 
quantity 5T{eV) ( |42[ ) is defined by exactly the same time 
integral as the interaction correction to the I — V curve, 
cf. Eq. (27) of Rcf. |^. In particular, in the limit of zero 
frequency and voltage one finds 



2L (1 _ e -t/floC)] is smallj with SS^o = 2Te 2 (3 J ta(t)e~ F ^ (l - e -^ R « c j dt. (44) 
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In accordance with FDT the combination in the right 
hand side is just the interaction correction to the zero- 
bias conductance of a coherent scatterertl multiplied by 
2T. 

We will now derive the interaction correction to the 
current noise in several important limits. 



D. High temperatures 

In the limit T ^> \/RqC it is sufficient to evaluate the 
function F(t) only at short times t < 1/T. In this limit 
from Eq. we get 



F(t) 



eH 2 
2C 



R 2R 2T 



(45) 



Expanding e in Eq. ( p2] ) to the first order in F we 
obtain 



SS U = 



e 2 /3 



(3R (eV eV 

1 + ^T (W^W- 1 



where function f(x) reads 

x cosh x 



2RCkQ 

uj±eV 



- 1 



2T 



(46) 



/(*) 



2 sinh x 2 sinh 2 a; 



(47) 



In the limit of small frequencies and voltages we then find 



e 2 f3 
ZRCh, 



(48) 



At high frequencies to ^> T, eV or large voltages eV ^> 
T, uj we get 



e 2 (3 



2RCkVL 



(49) 



Both results ( J4q ) and (B9) describe partial suppression 
of the current noise by Coulomb interaction. As we havg 
already discussed, Eq. ([48|) is consistent with the resultstl 



combined with FDT, whereas Eq. (49) just corresponds 
to the Coulomb offset AV^ = -e/3/2C on the I — V 
curve of a coherent scatterer at large voltages. For the 
sake of completeness we also note that in a specific limit 
\uj ± eV\ <C T Eq. ( [46| ) yields positive correction to the 
current noise 



6S„ = 



e 2 p 2 R \eV\ 
2ACTR 2 kVL ' 



(50) 



However, the magnitude of this correction is small in the 
parameter ~ e 2 R /(RCT) < l/g . 



E. Low temperatures 

Now let us consider the limit of low temperatures T <C 
\/CRq. At low voltages eV <C 1/RqC and times much 
longer than 1/R C the function F(t) reads 



^ 2 / smh[-KTt] 
W ~ ^ n [smhlnTRoC] 

Combining this expression with Eq. (|3 
small w, eT^ < T we obtain 



5 W =2T 



1 
Rk 



-(l-(TRoC)-> 



(51) 

in the limit of 

(52) 



For T > goEc exp(— go/2) the result ( p2\ j can be ex- 
panded in 2/g . In this limit for the interaction correc- 
tion we get 



8S*. 



In ■ 



1 



4T/3 

goRn RqCT 



(53) 



In the opposite limit of very low T < goEc exp(— go/2) 
(but still T ^> u, eV) the last term in ( |52] ) can be ne- 
glected and the interaction correction becomes 



2T/3 
~Rk~' 



In the limit T«w,ey« 1/RqC we can set T 
Eq. @. Then we obtain 



(54) 
in 



2Rk ^ 



1 



Rk 



R + Rs 
•±eV\[\oj±eV\R C] 



(55) 



If both u and V tend to zero, the last term in (|35|) can 
again be neglected and we find 



6S„ = 



HP 

Rk ' 



(56) 



If, however, uj and/or eV exceed the scale 
,g £ ; cexp(-go/2), one expands Eq. (J^) in 2/g 
and gets 



In - 



1 



2M/3 

5o^?k |w|-RoC 



(57) 



This expression applies for u> 3> eV. In the opposite limit 
in Eq. ( |57| ) one should simply substitute instead of uj. 
Note, however, that in the latter limit the corresponding 
result yields only one contribution (SS^) to the inter- 
action correction. Another contribution (SS^) will be 
found in Sec. 3F. 

To complete this subsection let us find the interaction 
correction in the limit uj,eV 3> 1/RqC. At large volt- 
ages the dependence of F(t) on V should be taken into 
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account. Evaluating the corresponding (linear in V) cor- 
rection to Fit) (fiT|), we obtain 



(58) 



These results complete our analysis of current fluctua- 
tions in coherent conductors with electron-electron inter- 
actions. 

IV. SUMMARY 



where 

w(x) = — 2 — 2x arctanx — In a; 2 + ln(l + x 2 ) + \x\tt. (59) 

For \ui ± eV\ ^> 1/RqC the asymptotics w(x 3> 1) ~ 
1 jZx 1 should be used. In this case we again recover Eq. 
(p9|). If, however, |w±eV| < 1/RqC, then the interaction 
correction is governed by another asymptotics w(x <C 
1) ~ —2 — In a; 2 and, hence, this correction is positive. 
Such an increase of the noise at \u> ± eV\ < 1/RqC is 
similar to that found at higher temperatures. 



F. Large voltages 

Now let us evaluate the remaining correction 8S^ 
(pq). At high temperatures T ^s> 1/CRq we obtain 



5S™ 



e 2 7 / eV 



e27 -0(|ey|-M)tanh |ey| - M 



RCnfl 



2T 

if w, eV > T. 



(60) 
(61) 



In the limit w, eV < T < 1/Ci? one finds 



= 



In- 



1 



2 7 (eF) 2 

3Tg RK *" Ti? C 



(62) 



Finally, at higher frequencies and voltages w, ey >Twe 
derive 



55 



(7) 



<7o RkQ 



+^tMln[l + ((|ey|-HW 



(63) 



Note that the correction SS^ is positive in all cases. 
As compared to previously obtained contribution SS^ 
the correction (63) becomes important in the limit eV 3> 



T, u. For such voltages both corrections add up, SS^ — 
SsL m +SS { J\ and yield 



_2C8-^QjeVl ln 



1 



g Rn ~ \eV\R C 
ifT, \u\ < \eV\ < 1/Ci? , 

SS " = - ^^T £ > if|e^|>T,l/CJ2o,M 



, (64) 



Combining the standard scattering matrix approach 
with the effective action formalism we have analyzed the 
effect of electron-electron interactions on current noise in 
mesoscopic coherent conductors in the metallic limit (]?]) . 

We have found that Coulomb interaction always leads 
to partial suppression of the Nyquist noise. The corre- 
sponding interaction term is proportional to the param- 
eter (3 (|l|) similarly to the interaction correction to the 
conductance^. Interaction-induced suppression of both 
conductance and Nyquist noise has the same physical ori- 
gin, and a direct relation between these two effects can 
easily be established with the aid of FDT. 

The effect of electron-electron interactions on the shot 
noise is somewhat more complicated. In this case we 
have recovered two different interaction corrections en- 
tering with opposite signs. One of them is negative and 
it is again governed by the parameter (3. Another correc- 
tion is positive and it is proportional to the parameter 7 
(^) which depends on the transmission distribution in a 
different way. The net interaction correction to the shot 
noise scales as 8S u {y) cx 27 — (3, i.e. it can be both neg- 
ative and positive depending on the relation between j3 
and 7. The contribution to 8S U (V) from the n-th con- 
ducting mode is negative provided its transmission T„ is 
smaller than 1 /2 and it is positive otherwise. For coherent 
diffusive conductors 27 — f3 — —1/15, i.e. in this partic- 
ular case electron-electron interactions tend to decrease 
the shot noise. 

The presence of two interaction corrections to the shot 
noise has a transparent physical interpretation. The (3- 
correction is due to Coulomb blockade suppression of the 
Landauer conductance (|^) while the 7-correction origi- 
nates from, |tbe term — T 2 in the expression for the 
shot noiseBlm The absolute value of this term is also 
decreased by interactions. But, since it enters with the 
negative sign, the corresponding contribution to the noise 
spectrum turns out to be positive. We believe that the 
effect of electron-electron interactions on higher cumu- 
lants of the current operator can be described in a similar 
manner. 

This work is part of the Kompetenznetz " Funktioncllc 
Nanostructuren" supported by the Landestiftung Baden- 
Wurttemberg gGmbH. One of us (A.V.G.) acknowledges 
support from the Alexander von Humboldt Stiftung. 
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APPENDIX A: EFFECTIVE ACTION AND OBSERVABLES 



Following^ let us combine the effective action formalismll§ll3 with the usual Landauer scattering approach. Within 
the latter approach one introduces a (relatively small) scatterer which connects two bulk reservoirs. The scatterer is 
described by the scattering matrix. In order to include electron-electron interactions it is necessary to reckon with 
the many-body Hamiltonian 



H = / dr* + (r) 



2m 



W(r) 



*(r) 



- / dr ( dr'^+(r)^ + (r')—^- 



-tf(r')*(r) 



(Al) 



Here the term W(r) accounts for boundary and impurity potentials, external fields etc. After the standard Hubbard- 
Stratonovich decoupling of the interaction term one arrives at the following path integral over an auxiliary field V(r, t) 



-iHt 



fVV(r,t>) 



•f e -iJ*dt'H e{t [V(.r,t')] 



J VV(r, t')e*/o dt 'I dr — i 

Here T is the time-ordering operator and 

H eff [V{r,t')] = [ drV+(r 



(A2) 



2m 



+ W(r)-eV(r,t) 



*(r). 



(A3) 



We choose to define the electron charge as — e. 

The time dynamics of the density matrix p is described by means of the evolution operator J defined on the Keldysh 
contour. In what follows we shall denote the field V on the upper and lower parts of this contour by Vi -2 . The general 
expression for the density matrix reads 



p(Xif,X2f,tf) — J dX\idX%iJ (Xif,X2f,Xu, X 2i ; tf, ti) p(Xu,X2i, U), 



(A4) 



where X the set of relevant quantum degrees of freedom. We shall assume that interacti on with the fluctuating fields 
Vx t 2 is turned on at a time i, = 0. The time tf is supposed to be large. Making use of (A2) and integrating over the 
fermionic degrees of freedom, we obtain 



J = J VViDVzexpiSlVuVi], 



where S is the effective action 



iS[Vi, V 2 ] = 2Tr lnG^ 1 + i% f dt'M 

* Jo 



LRl 



V, 



LR2\i 



(A5) 



(A6) 



Here we defined Vlrj = Vhj — Vrj and neglected the spatial dependence of the fields Vli,2 and Vm,2 inside both the 
left (L) and the right (R) reservoirs. The Green-Keldysh matrix Gy(Xi, X 2 ) (here X = (r, t)) obeys the 2x2 matrix 
equation 



iA.i-H (r 1 )i + eV(X 1 ) j G\ ( \\ . X . ) = ,m .V, - X,)a.. 



(A7) 



where Hq — (— V 2 /2m) + W(r), V is a diagonal 2x2 matrix with components Vij = ViSij and a z is the Pauli matrix. 

The above equation for the Green-Keldysh function should be supplemented by the initial condition for the density 
matrix p(t = 0) = Po, where p is the equilibrium density matrix of noninteracting electrons. In what follows we will 
need the solution of Eq. ( A7) in the absence of the fulcutating fields V1.2, which reads 

Gu(h,t 2 ) = -^0(^l-^2)^7l(<l,^2) + ^^/l(<l,o) / 5o^7l(o,^ 2 ), 

G 22 {t x M) = -i6{t2-t 1 )U 2 {t 1 M) + iMti^)poU 2 ^M), 
G 12 (h,t 2 ) = iUi(t u 0)p Q U 2 (0,t 2 ), 

G 21 {t u t 2 ) = -iU 2 (h,0)li~p Q ]Ui(0,t 2 ), (A8) 
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where U\^ 2 (ti, t%) are the evolution operators 

Uifi{h>h) = T exp 



i J df (H -eV h2 (r,t') 



(A9) 



One should keep in mind that in the operator products like UpU integration over intermediate spatial coordinates is 
implied. 

Instead of specifying W(r) we will describe electron transfer between the reservoirs by means of the scattering matrix 
formalism. This procedure is standard and further details are provided in Appendix B. In calculating the trace in 
Eq. ( |A6| ) , we shall make an explicit integration over the longitudinal coordinates. Integration over the transverse 
coordinates is replaced by summing over the transmission channels of the scatterer. It is convenient to introduce the 
phase variables 



<p+(t)= / dt'(eV LRl (?)+eV LR2 (?))/2, 



<p~(t) = / dt' (eV LR1 (t') - eV Lm (t')) 
Jo 



(A10) 
(All) 



Provided the dimensionless conductance go is large j > 1 fluctuations of f~(t) are strongly suppressed, so we can 
expand the exact action zSJ^] in powers of ip~ keeping the full nonlincarity in <p + . Note, that the external voltage 
enters only in ip + but not in tp~ . Hence, for the system of Fig. 1 we have to impose the constraints Ylj Vjit) = e Kr^> 
fj{t) = 0- Here and below in this Appendix the summation runs over different elements in our circuit, i.e. the 
scatterer and the shunt. 

Let us define the kernel of the current operator. Consider, e.g., the upper part of the Keldysh contour. The 
charge is proportional to the derivative of the action with the respect to the external voltage. Due to the property 
SS/5Vx(r,i) = e J 5S/Sipi(r,t) we can relate the difference of the charges of the left and right reservoirs to the 
integral of 5S/SipL R (t) over time. Hence, the latter derivative is proportional to the current through the scatterer. 
Keeping track of proper ordering in the upper part of the Keldysh contour, for t > ? we obtain 



(i 3 {t)i 3 {t')) = -e 2 / V<p 



8ipij(t)5tpij(?) 



Combining (A12) with a similar relation in the lower part of the Keldysh contour we get 



Sipij (i)6ipij (?) 8ip 2 j (t)6ip 2 j (?) 



(A12) 



(A13) 



This equation is equivalent to (|1^ 

Below we shall proceed with an explicit calculation of the action by defining the scattering states. 



APPENDIX B: TRANSMISSION CHANNELS 

Following the usual S'-matrix approach^ let us introduce the transmission channels. We will assume that far from 
the scatterer the electron propagation in transverse and longitudinal directions can be described separately. In this 
case the solution of the Schrodinger equation 

V 2 

t/;(r) + W(r)iP(r) = Eip(r) (Bl) 

2m 

can be factorized 

^(r) = Cn$n(r±)Xn(x). (B2) 

Here x is the coordinate along the lead and r± are the transverse coordinates. The transverse wave functions <5 n (r_i_) 
satisfy the equation 

V 2 

-77 ± ^n(r±) + W(r ± )$ n (r ± ) - E n $(r ± ), (B3) 



11 




'L "R 

FIG. 2: Scattering states 



where the subscript n enumerates the transmission channels (we are considering only channels with E n < Ef)- The 
function Xn{x) outside the scatterer region is defined from the equation 



2m dx 2 



Xn{x) =(E- E n )Xn{x). 



(B4) 



Since the electronic states with energies E close to the Fermi energy Ep mainly contribute, it is sufficient to describe 
the elec tron dynamics q uasiclassically. We define the energy £ = E — Ef and the particle velocity in the n-th channel 
v n = y/2(Ep — E n )/m. Then the wave function can be expressed as 



Xn{x) 
Xn{x) 



—imv n x jpout 



f° ut {x), left reservoir. 



n imv n x „out 



° ut (x)+e tmVnX g™ (x) , right reservoir. 



(B5) 



In this way we have introduced the envelopes of the fast oscillating functions cxp(±imv n x) . Consider first the left 
reservoir. The functions f™(x) and f° ut (x) satisfy the following quasiclassical equations 



-iv n —f™{x)=ZC{x), 
iv n -^f™ t (x)=Zf™ t (x) 



with the solutions 

Analogously, for the right reservoir we find 



i£ x I v n 



c / v n 



p —i£x/Vm p i£x/v n 



I V 7n -yj Vm 

The eigenfunction of the whole system with the energy £ in the left reservoir may be expressed as 

^( r ) = E [aLne imVnX ft(x) +b Ln e- imv " x f™ t (x)] $„(r ± ), 

n 

while in the right reservoir we get 

M r ) = E [ b Rk^ mVkX 9k Ut (x)+a Rk e- m ^gT(x)] * k {r ± ). 



(B6) 



(B7) 



(B8) 



(B9) 



(BIO) 



The amplitudes of the outgoing r and incoming r states (see the Fig. 2) are related via the scattering matrix 
5(0 



( bhi \ ( am \ 



^LN L 

b R i 



V bRN R I 



0,LN L 

am 



(Bll) 



V aRN R ) 
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The unitary matrix S with dimensions (Nl + Nr) x (Nl + Nr) has the block structure 

*o-(M)" 

The diagonal blocks f and f' describe reflection back to the left and right reservoirs, respectively. The off-diagonal 
blocks describe transmission through the scatterer. Later we shall neglect the ^-dependence of S. 

Let us now combine the incident f™(x) and outgoing fn Ut {x) wave functions belonging to the same channel into 
one wave function i/j n (x). Namely, we assume that the scatterer is located at x — 0, and for the left reservoir [x < 0) 
we put 

^-{P(-,),',to. < B13 > 

Analogously for the right reservoir (x > 0) we define 

These new functions are defined in all the range y £ [— oo, +oo] and are equal to 

^(y) = (B15) 



Let us emphasize, that here the index j enumerates all conducting channels, both in the left and in the right reservoirs 
(for convenience, we assume that the left channels are enumerated first). 

In the presence of the fluctuating field V(t), the matrix elements of the Hamiltonian in the reservoirs take the form 

d 

Hij = -iviSij — - eVi(t)6ij, (B16) 

where Vi — Vl for all left channels and Vi — Vr for the right channels. If at initial time t\ the wave function was 
i/> n (ti,y), then at the final time t 2 > t\ it becomes 

M*2,V) = e^W-v^Hntyuy-Vnfo-tx)), y<0 or y > v n (t 2 - h); (B17) 
Mh ,y) = V e *™^)-^^)]-K(^-^)-^(*-^)] Sn)fc & fe U, % - v k (t 2 - h)) ,0<y< v n (t 2 - tl ). 

Here s nk are the matrix elements of the 5-matrix and, as before, <p n (t) = L dieV n (t). On the other hand, by definition 
of the evolution operator we have 



i>n(t 2 ,y 2 ) = / dyiU n k{t2,h;V2,yi)i/>k{ti,Vi)- 

7„ •* 



Comparing this expression with (B17), we find 

U n k 



(t 2 ,h-y 2 , yi ) = e *»w\^5 (V*—^ -b+h) + 6(y 2 )9(v n (t 2 - h) - y 2 )e^"( t2 -^) [s nk - 5 nk ] x 



It is convenient to introduce the new coordinates r = y/v n . More precisely, instead of the wave function with 
the components ip n {y) we introduce the functions rj n {r) = y/v^ip n (y/vn)- The kernels of the operators will also be 
transformed. If the two functions are related to each other by means of a linear operator 

^ 2) (2/) = E/ d 2/ KnkM^iv'), (B18) 
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then the corresponding wave functions r]^ and t/ 1 ' satisfy the following relation: 

*i 2) M=£ j dr' k nk {r,r , )^\r'), 

k J 

where 

K nk {r,r') = ^/v^K nk (v n T,v k T')y/vk~. 
In this representation the evolution operator can be simplified. We find 



(B19) 



(B20) 



u(t2,h-,T 2 , n) = s(t 2 -r 1 -t 2 + hy^U i + e^ei-r^e-^-^is - ly^-^ L-^l (B21) 



The matrix (p is diagonal with respect to the channel indices (pi k = ifiSik- We also obtain an expression for the inverse 
operator, i.e. the operator defined by J dr 2 U (t 2 ti; t^t^U^ 1 (t 2 ti; t 2 ti) = 5(t3 — ti). It reads 

^ 1 (i 2 ,ti;r 2 ,r 1 ) = S( n -T 2 -t 2 +t 1 )e^{i + 6(r 1 )9(-T 2 )e-^- T ^[S + - i] e ^- T Ae~^ . (B22) 



Eqs. (B21), ( B22 ) apply for t 2 > <i, in order to construct the corresponding expressions in the opposite case one 
should just use the property U(t 2 ,ti;r 2 ,Ti) = U~ 1 (ti,t 2 ; r 2 , n). 

Finally let us define the equilibrium density matrix for noninteracting electrons. It can be written in the form 



Po,nk{yi, 2/2) = 8 nk 



dp e 



ip{vi-V2) 



2tt l + eP v -/ T 2 



5nk s , s Snk ^iT 

<%1 - V2) 



2tt 



v n sinh 



Performing the transformation (B20) we obtain 



po{n,r 2 ) = - 5(n - t 2 ) 



iT 



sinh \kT{t\ — t 2 )] 



(B23) 



(B24) 



APPENDIX C: EXPANSION IN THE PHASE DIFFERENCE 



We shall expand the effective action (A6) perturbatively in ip~ . The field ip + will be taken into account exactly in 
each term of this expansion. The expansion starts from the first order in tp~ , since for (p~ = the contributions from 
the forward and backward parts of the Keldysh contour cancel each other. We get from Eq. ( A7) 



2Tr InG 



v 



2Tr In 1 



G (p- 



(Cl) 



The Green-Keldysh matrix Go is evaluated for ip" — 0, i .e. it is defined by Eqs. (A8) with the evolution operator 
(B21) taken at (p — (p + ). The fluctuating field (p~ in (CI) is a unity matrix in Keldysh space and a diagonal matrix 
in the space of conducting channels. Performing the summation over the Keldysh indices we arrive at the first order 
in ip~ contribution to the action iS^ 



iS^i^] =-i dt I ds I dn / dr 2 Tr 



, 1 iriT 



U- 1 (t,0;T2,s)^-(t) 



(C2) 



7T sinh[7rT(ri — t 2 )] 

For simplicity in Sec. 3 we have set tf — > 00. Here we will keep it finite and use the conditions <^~(0) = <p~(tf) = 0. 



The 5-functions contained in ?7-matrices of (C2) will lead to a singularity of the form 1/ sinh[7rT(ri — t 2 )) which is 
cured as follows. Let us choose close but not exactly equal arguments s± and s 2 . Expanding the combination 



dt / dr\ I dr 2 Tr 



U(t,0;s 1 ,T 1 )U- 1 (t,0;T 2 ,s 2 )<p-(t) 



to the first order in si — s 2 and multiplying the result by 1/ sinh[7rT(si — s 2 )], we obtain 



(C3) 



(C4) 
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Making use of the condition Tr = Tr we get 

iS {1) [^] = -~Tr[i + i\ J d T ip-(T)<p + (T), 



where if ± {r) = (p±(r) - <Pr(t). 

Consider now the contribution to the action of the second order in ip~ . It is defined 



(C5) 



as 



iSWfo*] = -Tr 

After a straightforward algebra one obtains 

*/ t, 

i5 ,(2) [< / 3 ± ] =- / dri J dT 2 p (T 2 - Ti)po(n - t 2 ) : 



Gi2<p G 2 np~ 



(C6) 



o o 

Tr j e -<0 + (*i)-$ + fa)] [S+0-( n )S - ^-(T 1 )]e l[ * +(ri) " ¥5+(r2)1 [S+£-{t 2 )S - <^(r 2 )] \. (C7) 
Taking into account the block structure of the S'-matrix, we find 



iS^l^} = -±-J d Tl J dr 2 a{n -r 2 )i [Tr (t' + f) 2 + Tr (t+t) 2 ] <p~ (n)<P~ (n) + 

n n 



2Tr 



o o 
f'r' + ii+ 



; (<p + (n) - ip + (r 2 )) ip (n)<p (r 2 ) 



(C8) 



where a(r) is defined in Eq. (pT|). Introducing the parameter /3 = Tr 
more compact form 

t, t, 

Tr \t+t 



r'r ' tt+ 



/Tr we rewrite Eq. (p| in a 



■.S^l^} = Jdn I ' dT 2 a{Ti-T 2 )y-{Ti)y-{T 2 ) [l-/3 + /3cos (^ + (n) - V + (t 2 ))] 



(C9) 



o o 



We now proceed to the third order contribution iS^ to the effective action. It reads 



iS^ = iTr 
2 



dti j dt 2 j dt 3 6(h - t 2 )8(t 3 - t 2 )T{h)T{t 2 )T{h)p^ 

.0 

Here we used the notations 

iT 



i-Tr 
12 



(Gi 2 <fi + G 2 \(p~ 



P?(n - r 2 ) = - - sinh[7rT(ri _ T2)] i, A*) - U-\t,0)t-W(t,0)- 



We obtain two terms from Eq. (CIO). The first one 



iSf^} = ^Tr[£+f] / dr (^-(r)) 3 ^+(r) 



67T 



(CIO) 



(Cll) 



(C12) 



u 

In deriving this result we employed the same - although somewhat more involved - regularization p roced ure as for 
the first order contribution This procedure allows to determine the correct overall prefactor in (C12). One can 
then verify that the resulting effective actio n satisfies the requirements of FDT. 
The second term, coming from Eq. (C1C), has the form 



-Tr 



\ii + ) 2 r'? + 



f-tf f-tf ftf 

I dyi / dy 2 / dy 3 p is {y 2 - yi)Po S (y3 - V2)Po S {yi ~ ys) X 

JQ JO 



" ' JU JU 

<P~ ~(yi)v~ ~(y2)<P~ "(2/3) {sin(ip + (y 2 ) - f + (yi)) + sin (</? + (y 3 ) - ip + {y 2 )) + sm((p+(y 1 ) - (p + (y 3 ))} ■ (C13) 
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Defining the parameter 7 = Tr 
now all four contributions ( |C5 



/Tr [it + ] and shifting the phase ip + by eV we obtain Eq. (j35|). Collecting 
we arrive at the final result for the effective action 



), flClg) and (|C1 

s = s {1) + s {2 



c(3) 
P 



s< 3 >. 



(C14) 



This action is valid up to the third order in ip , and the variable (p + is treated exactly in each of the terms in (C14). 
It is instructive to compare our results with the AES actionta derived for tunnel junctions (fi — > 1) to all orders in 



Rewriting the actiorJij in our notations together with the capacitive term one has 

tS AES = -j-J dtxj ^ 2 a / (t 1 -t 2 )0(t 1 -< 2 )sin(^+(t 1 )-^ + (t 2 ))sin^-^cos-^^ 

dt 2 a{ti - t 2 ) sin V ^ sin V ^ cos((^+(ti) - V + (t 2 )) - 4 / ' dtC(p 
2 2 e. h 



f ft ' 
dti 
Jo 

2 f tf f tf 

5p/. *W 



(C15) 



e 2 i? 



Here we denoted aj(f 1 — t 2 ) = <5'(ii — £ 2 ). This <5-function should be understood as a smeared one. 

Let us expand (C15) in (p~ and compare with our results order by order. The first order terms are exactly the 
same for both models. The difference between the models shows up in the second order terms, for our model the 
parameter (3 appears exp licitly in the second order cont ribution ( p9| ). In the limit (3 — > 1 this expression reduces to 
that obtained from (C15). Expanding the action (C15) to the third order in ip~ one only recovers the term of the 
form QC12 ) with (3=1, while another term (C13) cannot be recovered. Contributions of this nature are not contained 



in the AES action at all since they are proportional to higher orders of the channel transmission T n , 

It is worthwhile to point out that a formally exact pec*£seHtation for the effective action of a coherent scatterer (all 
orders in T n and all orders in <p ) can also be derivealj'Ej'Prr 1 !. However, this formal expression turns out to,^ quite 



complicated to deal with in the situation addressed here. For g = R q / R ^> 1 and provided instanton effectsEilO can 
be neglected all necessary information is equally contained in a much simpler form of the effective action derived in 
the present paper. 

We also note that there exists a simple relation between the action derived here and the cumulant generating 
function describing the full counting statistics of the charge transport in noninteracting coherent conductors^. This 
relation can be established if one neglects fluctuations of the phase variable, i.e. sets tp + = eVt, and chooses ip~ to 
be time independent. By identifying tp~ = —A and expanding the generating function ln(x(A)) (defined in Eq. (37) 
of Ref. p3) one arrives at the following indentity: 



ln( X (X)) = -S[eVt,-X}. 



Finally, let us use the AES action (C15) in order to illustrate the imp ortan ce of the third order in ip~ terms for the 
calculation of the currect-current correlation functions. App lying Eq. ( A13 ) one gets a contribution stemming from 
the double differentiation of the term with ai in the action (C15) 



.tfl / . /' ) = - ( - t') ai (t - t') sin(/ (t) - v+ (£')) sin y ^ - 9 ^ ) + ( / - - t'). 



(C16) 



Evaluating this average with the aid of path integrals one has to keep all nonlinear terms in the pre-exponent. However, 
the dependence of the cosine term on ip + in the action in the exponent can be neglected provided j > 1. Applying 
the identity sm(ip+(t) - tp + {t')) = X) v= ± e w ^ +(t)_ ' <9+(t ' )1 /2j we arrive at the following integral 



XV exp <iu(ip + (t) -ip + {t')) 



dtcp + (t) 







C(p~{t) 



R 



which yields the (5-function 



2tt 



ip~(t) = — v 6(t-t) 1- e 
9 



-{t-^/RC^ _ Q{t , _ ~ t) ^ _ e -(t'-t)/RC 
I 



(C17) 



(C18) 



As it was expected, <p (t) is indeed small for g ^> 1. 
Combining the above expressions we arrive at the contri- 



bution 



SS{t, t') 



2 5{t-t')/2RC. 
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The same contribution multiplied by the factor /? was 
derived in Sec. 3B from the term Sp (|3l"l). The above 
analysis makes the significance of the third order in ip~ 



terms in the action particularly transparent: The ker- 
nel aj{t) introduces the time-derivative of e - t / RC which 
compensates for an extra smallness ~ 1/g. 
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